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I. INTRODUCTION 


In gauge theories at a nonzero temperature T, the order parameter for deconfinement in 
a SU(N C ) gauge theory is the Polyakov loop, 


L(x) 


-l/T 

-^rtr'Pexp \ig j A°(x,t)cIt 


( 1 ) 


This loop is related to the propagator of an infinitely massive test quark, which sits at a 
spatial point x, and just propagates forward in imaginary time, r : 0 —> l/T. A quantity 
invariant under small gauge transformations is obtained by tracing over color, denoted by 
tr; V represents path ordering in r, g is the gauge coupling, and A 0 the vector potential in 
the r direction. 

To be definite, we consider the case of a test particle, and its associated loop, in the 
fundamental representation of a SU ( N c ) gauge theory. (As for test particles, though, there 
are (independent) loops for all irreducible representations.) For two colors the loop is real, 
but for three or more colors the anti-loop, L*(x), is an independent quantity, representing 
the propagation of an anti-particle backwards in (imaginary) time. 

The correlation functions of Polyakov loops are fundamental to our understanding of the 
behavior of gauge theories in a thermal bath. These have been studied in great detail, both 
perturbatively, and from numerical simulations on the lattice [1-5]. In general, understand¬ 
ing the behavior of the loop correlation functions is involved, requiring the use of effective 
Lagrangians and non-perturbative input, such as from numerical simulations on the lattice 
[ 6 ], 

This is not true in an Abelian gauge theory, where at least at tree level, one can imme¬ 
diately compute the correlation function between two loops. In this case, the path need not 
be straight but can be arbitrary. A source couples to the photon as ~ exp(f <f ds J^A^), 
where “s” the worldline of the test particle. It is then trivial to integral over over A 

ds cf ds' r 




and obtain the sources tied together by a photon propagator, A fiv . 

This exercise is not possible in a non-Abelian gauge theory, where the fields are matrices 
and do not commute. As an exercise, in this paper we compute the correlation functions 
between Polyakov loops at tree level. Even in this instance the generalization of Eq. (2) 
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is not possible for arbitrary N c . This is basically because for 3 < N c < oo, one cannot 
determine the element of the Lie group, the Wilson line, from that of the Lie algebra, A 0 . 
We compute the two point functions of loops in two special limits, for two colors, and in the 
planar limit, N c —> oo with g 2 N c held fixed. 

As a computation at tree level, our results are valid only in the limit when the temperature 
T is very large, so that the running coupling constant, g 2 (T), is small. We compute the 
correlation functions between two loops at an arbitrary distance, x; this will always appears 
with the only dimensional scale in the problem, the temperature T, as xT. At large distances, 
this must agree with the results of an ordinary perturbative expansion computed to tree level. 
Because each loop is a trace over color indices, the leading diagram, from one loop gluon 
exchange, cancels, and so unlike the Abelian theory, there is no term ~ g 2 / ( xT ). The leading 
diagram is given by two gluon exchange, which is rsj g A /(xT) 2 . This is valid at both N c = 2 
and N c = oo, and is a trivial check on our results. 

The natural expansion parameter in the correlation functions is ~ g 2 /(xT). Thus at 
short distances, naive perturbation theory cannot be used. When the loops are very close 
together, though, xT <C 1, though, physical intuition suggests that it shouldn’t matter how 
the loops are tied off from one another. Thus the result should be the same as for a Wilson 
loop that is very long. As a Wilson loop is a single trace over color, single gluon exchange 
contributes, so the result at small distances should be Coulombic, ~ g 2 / (xT). 

We find that at short distances, in both cases the leading behavior at short distances is 
the same for the correlation between a Polyakov loop and an anti-loop (or two Polyakov 
loops), and the corresponding Wilson loop. For two colors, at short distances the behavior 
is Coulomb-like, although the coefficient is not identical to that from single gluon exchange. 
For an infinite number of colors, the behavior is milder, ~ \/g 2 / (xT). Although our results 
are only valid to tree level, for arbitrary distances we obtain non-trivial functions of g 2 /(xT) 
which can be written in closed form. 

In the Conclusions, Sec. (V), we discuss why the non-Coulombic behavior at short 
distances is an artifact of the planar limit, and discuss how to compute the short distance 
behavior from effective theories [3, 4], Our study was inspired by recent studies of the 
susceptibilities of Polyakov loops by Lo, Friman, Kaczmarek, Redlich, and Sasaki [5]. We 
also discuss in the Conclusions the implications of our analysis for the measurements of such 
loop susceptibilities. 
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II. CONVENTIONS 


Since we only compute Polyakov and Wilson loops, at the outset it is convenient to 
adopt static gauge, <9 0 A°(t, x) = 0. It is also necessary to fix the gauge dependence of 


spatial gluons, A\ at a given time, but at tree level the spatial gluons do not enter into the 


correlators of Polyakov loops or of the Wilson loops which we consider. 

The advantage of static gauge is that for the Polyakov loop, time ordering can be ignored, 
with the Polyakov loop just the exponential of a single field A°(x). For single gluon exchange, 
it is clear that in an arbitrary gauge, time ordering gives the same result as in static gauge. 
It is not so obvious for the exchange of two or more gluons, but in the end, both Polyakov 
(and Wilson) loops are gauge invariant. In any case, what we are really interested in is how 
the correlator of the exponential of non-Abelian A 0 fields is related to the correlator of the 
A 0 fields. 

The gluon field A 0 = A° a t a , normalizing the generators as tr (t a t b ) = 5 ab /2, a, b = 1... N% — 
1. For any N c , in static gauge the gluon propagator for A G a is color diagonal, 



( 3 ) 


Since we only compute two point functions, which are a function of a single spatial distance 
x, henceforth we replace x by x = |T|. We introduce A (x)/T since Polyakov loops are 
dimensionless, and so at tree level their correlation functions are functions only of xT. 

As our computation is only at tree level, several well known physical effects do not enter. 
Beyond leading order, the Debye mass, ~ gT, arises to screen correlation functions involving 
A 0 . Also, the renormalization mass scale enters to ensure the coupling constant runs. 

The virtue of computing at tree level is that it is mathematically well defined. Of course 
then one must be careful to consider the limitations of this approximation, as we discuss in 
the Conclusions, Sec. (V). 

III. TWO COLORS 

For two colors, we can use the well known property of Pauli matrices, cr a , and evaluate 
any element of the Lie group from that of the Lie algebra, 


exp (i (3 a u a ) = cos (/3) +i (3 a cr a sin (/3) , 


( 4 ) 
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where = p a /3, 0f = 1. 

The vector potential is A 0 = A®cr a /2. Since A° are the only degrees of freedom which 
enter, we denote it simply as A a , and use a convention for repeated indices in color, 
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Al = E(4,) 2 . (5) 

a= 1 


Hence the Polyakov loop is given by the first term, from the cosine: 


L{x ) = - tr 


JgA° a (x) t r a /(2T) 


= COS 


oo 

= E 


(-l) n ( 9 


^ (2n)\ V4T 2 

n =0 v ’ v 


A a (x)' 


( 6 ) 


This also shows that for two colors the loop is always real, L = L*. The correlation function 
of two loops is then 


(L(x)L( 0)) = ( E 


(- 1 )" ( 3 


~Aa{x) 




( 7 ) 


^ (2 n)\ \4T 2 w / ^ (2n')! \4T 2 

We ignore self energy corrections, which vanish with dimensional regularization. Then the 
correlator only receives contributions from terms where n = n', 


(mm) = e 


n=0 


((2n)!) 2 V2 T) 


%T <HW 2 )” H(0)T> 


( 8 ) 


If only one color contributed, then this correlation function would follow immediately. 
Because for two colors there are three gluons, we have to solve the problem of how many 
ways to tie three fields together in all possible ways. We do this in two ways, as a cross 
check on each method. The first method is combinatoric, just a matter of counting different 
fields together. The second is analytic, and may be useful for other problems, when counting 
permutations might be confusing. 


A. Combinatoric analysis 

We solve a more general problem, that of computing the combinatoric factor in Eq. (8) 
when the sum over fields is from a = 1 to m, instead of m = 3. We work inductively, 
starting from m — 1, where the answer is immediate, and then generalize to arbitrary m by 
induction. 

The problem is to determine the number of ways in which we can tie the fields 

<(A Q (a;) 2 )" (A(0)T> (9) 
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together. We denote this as D(2n,m). 

For a single component field, m — 1, the answer is obvious. The first A(x) can connect 
with any of 2n v4(0)’s; the second A(x) can connect to any of the remaining 2 n — 1 A(0)’s, 
and so on. This gives the usual factorial, T>{2n, 1) = (: 2n)(2n — 1)... = (2n)!. 

To solve for m > 1, we start with some special cases. 

When n — 1, the first A a {x) can connect to one of the two A fe (0)’s. The second A a (x) 
must then connect to the other A&(0). Since a = b, the sum over a gives m, and so 

T>(2, m) — 2\m. (10) 

Now consider n = 2. The first A a (x) can connect to any of the four Af,(0)’s. There are 
then two possibilities. The second A a (x) can connect to the partner of the first A*,(0), which 
gives a factor of m. Tying the second pair together gives J\f(l,m) = 2m, or 8m 2 in all. 

Alternately, the second A a (x) can tie to a different pair, A c (0) 2 . It can do so in one of 
two ways. This leaves two ways for the third Ad(x) to connect to one of the two remaining 
A(0)’s. All four A’s have the same index, so there is one sum over a, or 16m in all. The 
sum is 8m 2 + 16m, or 

T>{ 2,m) = 4!! (m + 2)m . (11) 

where n\\ = n(n — 2)... . 

For any n , the numbers of ways to connect with the highest powers of m is evident, just 
m n . The first A a (x) can connect with any of the 2 n Afe(0)’s. The second A a (x) has to 
connect with the pair of A&(0), which gives 2 n x m. The third A c (x) can connect to any of 
the remaining (2 n — 2) A^(0)’s, and so on. Thus for large m, the leading term is 

V(2n, m —» oo) = (2 n) m (2 n — 2) m.. .2 m — (2n)!! m n . (12) 


Looking at the results for n — 1 and 2, and the result for large m, a guess for arbitrary 
m and n is 

(2 n — 2 + m)H 


V(2 n,m) = (2n)!! 


(13) 


(m — 2)!! 

When m = 1 we define (—1)!! = 1, so V(2n, 1) = (2n)!!(2n — 1)!! = (2n)!. 

To establish the result for m > 1 we use induction, and consider how T>(2n + 2, m) is 
related to V{2n,m). As always, the first A a (x) can connect to 2 n + 2 different Aj,(0)’s. 
There are two ways in which the second A a (x) can connect to a field A^(0). The simplest is 
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if it connects to the pair of the A&(0) which the first A a (x) connected to. There is only one 
such Afe(O) to connect to. The sum over a = b gives a factor of m, so we have (2 n + 2 )m 
times the number of ways of tying n pairs of A a together, which is V(2n,m). 

Alternately, the second A a (x) can tie with a A c (0), which is in a different pair than the 
first A a (x) connected to. There are 2 n different A c (0)’s to choose from. This possibility 
looks complicated at first, but explicitly what we have is 

(A a (x)A b ( 0)) <A a (x)A c (0)} ((A c (*) 2 ) n A fe (0)A c (0) (^(O) 2 )"" 1 } . (14) 

Now A b (0) and A c (0) start out as being in different pairs. However, the propagators are 
diagonal in the color indices, and evaluating the propagators in Eq. (14) gives a = b and 
b = c. After doing so, we have that Ab(0) and A c (0) have the same indices. There is no 
factor of m, because we still have to tie Af,(0) and A c (0) to other A(x)’s. Consequently, 
we are left with (2n + 2) (2 n) times the number of ways of tying n pairs together, which is 
T>{2n, m). 

The sum of these two terms is 


V(2n + 2, m) = (2n + 2) (2 n + m) T>[2n, m) . 


(15) 


The solution to Eq. (10), for n — 2, Eq. (11), for n = 4, and this equation, is the general 
result of Eq. (13). 

For the case of interest, m — 3, we obtain an especially simple result, 


V(2n, 3) = (277.)!! (2 n + 1)!! = (2 n + 1)! . 


(16) 


With this symmetry factor in hand, from Eq. (8) the two point function of loops 


00 (9 I 1\| / 2 \ 2n 2i 

(L(x)L(0)) = Y, (( 2 W )!) 2 ' [h A j = ^h {z) + z sinh (z) , z = ^ — . (17) 


At large distances, z —> 0 or x —> oo, 
(L(x)L( 0)) ~l + ^ 2 + ^ 4 + •••-! + 
For general N c , to ~ g 4 [3, 4] 


3 (g 2 \ 2 i 5 
32 \4n) (.xT ) 2 + 24 


( g2 ) +... . (18) 
yl67ncT / v ; 


{L*(x)L( 0)) « 1 + 


N 2 C - l { g 2 \ 2 1 


8 N 2 \4n J ( xT ) 2 


+ ... , x —y oc . 


(19) 


7 








This result is due to the exchange of two gluons, with the factor of (A 2 — 1)/(2A 2 ) from the 
Casimir for the fundamental representation. This agrees with Eq. (18) for N c = 2. Notice 
that the term ~ z 4 ~ g 8 is of higher order than has been computed previously. This is 
because we only include tree diagrams. 

At short distances, z —» oo, 


1 


(L(x)L( 0)) « - exp(z) pa - exp + —-— + log 


9 


167 txT ) 


g- 


i-^0. 


16. xT ' - ' — ■ ■ ■ ~ * - < 2 °> 

The coefficient in the exponential is the same as for the Wilson loop at short distances, Eq. 
(41) in Sec. (Ill C). 

What is novel here is that Eq. (17) is valid for arbitrary distances, and interpolates 
smoothly between small and large values of g 2 /(xT). 


B. Analytic method 

We work with components, and use the multinomial theorem 

(A„(xf Y= V ”, Af'(x) Af'(x) Af°(x) 4.+ ta+t „„ . (21) 

fel,fe2,fc 3 =l Kl ' 2 ' 3 ’ 

After expanding each bracket in Eq. (8) we ignore self energy correction to match the powers 
of A 0 , 

o° i ^ 

< £ <*> i(0 >>=£«(£)’ (22) 

ki,k 2 ,k 3 =0 ' 1 2 3 ' 

Since we have evaluated the fields component by component, it is direct to evaluate the 
correlation functions, 

(A\ (x) 2kl A 2 (x) 2k2 A 3 (x) 2ks A 1 (0) 2fel A 2 (0) 2fe2 A 3 (0) 2fes ) 

= (Axixf^Axi 0) 2fcl ) (A 2 (x) 2k2 A 2 ( 0) 2fc2 > (A 3 (x) 2k ^ A 3 (0) 2fc3 > 

= (2 kx)\ (2 k 2 )\ (2 k 3 )\ A 2fcl A 2fca A 2fca = (2 kx)\ (2 k 2 )\ (2fc 3 )! A 2n . (23) 

The last equality was obtained by taking into account the relation k± + k 2 + k 3 = n. 

Next we represent the Kronecker symbol as 

X, , , — \ J) ^ /-k 1 +k 2 +k 3 -n /OA.\ 

0ki+k 2 +k 3 ,n 2 ^ JP ^ S > \^) 







where the contour C includes the origin. 
The two point function becomes 


1 °° I 

(L(x)L(O)} = — £ ^ )!)2 V4r 


" 2 A 


2 n 


d( 


c 




C " +1 \£W ! > 2 


(25) 


We then evaluate the sum 


E^fc‘ = 


+ C” + 1 e(C), 


k =0 




VT^4C 


(26) 


where F(() is holomorphic about the origin, ( = 0. Note the factor of £ n+1 in front of F((). 
Because of this, only the first term in Eq. (26) contributes, 


d( 


+ C n+1 F(0) =;r^ 


d( 


2 ni J c C n+1 (1 - 4C) m /2 • 


2vri Jc C n+1 V VI -4C 

Since the function F(Q is complicated, this helps to simplify the problem. For m = 3, 

1 _ 1 _ _L Y' 2 2 i + if ri 

C” +1 (1 - 4C) 3 ' 2 C” +1 jl r (I) ^ ■ 


(27) 


(28) 


Using the residue theorem, 


d( 


2vU Jc C n+1 (1-4C) 3 / 2 


= 2 


2 72 


r ("+!)) 

n!r (!) 


Altogether, 


Using 


(£M£(0)) = f; r( _ n l i) UA 


^((2n )!) 2 r(f) Vr 


2n 


r n+- = 


(2n + 2)! 
4 n+1 (n + 1)! 


7T , 


(29) 


(30) 


(31) 


it is direct to show that Eq. (30) agrees with the combinatoric result of Eq. (17). 


C. Wilson loop 

For completeness, in this section we compute a Wilson loop for two colors in the tree 
approximation. We take a Wilson loop which has the same sides in imaginary time, and 
just tie the ends at r = 0 and t — 1/T together. Since we drop self energy corrections, we 
can drop the terms from r = 0 and 1/T, and concentrate just on gluons exchanged between 
the two sides, at x and 0. 
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In general, there is no simple relation between the exponential of the sum of two elements 
of the Lie algebra, and the product of the corresponding exponentials. In this case, however, 
path ordering implies that they do factorize: 


V exp 



(A a (x)a a 



exp 



exp 



(32) 


There are several ways to establish this. One can check it to the first few orders in g, and then 
use an inductive proof in powers of A a (x) and Ab(0). A more elegant proof is the following. 
Let A a (x) be infinitesimally small, with Ab(0) arbitrary. Then by path ordering, A a (x) 
always lies to the left of any A&(0). For an inhntesimal A a (x), to linear order in A a (x) we 
have ~ (1 + ig/ ( 2T)A a (x)cr a ) exp(zg/ (2T)A b (0)a b ). We can then write a differential equation 
for A a (x), taking the derivative on the left. The solution of this differential equation is the 
product of the exponentials in Eq. (32). 

Thus the Wilson loop is given by 


W 



exp 



exp 


We can now use the relation of Eq. (4), 



A b (0)a 



(33) 


w= ( cos {4f^ Aa ( x ) 2 ) cos {4f^ Ab ^ 2 ')) 

+ ^ A a (x)A a ( 0) sin A 2 b {x)^ sin (°))^ • ( 34 ) 

For general N c , the Wilson loop factorizes as above, but then we cannot work out the 
expansion of the element of the Lie algebra. 

The Erst term in Eq. (34) is equal to the correlator of two Polyakov loops. The second 
term is new, but can be computed by similar means. In a power series expansion, this second 
term equals 

OO . —1—0 

£ ((2ra + l)!) 2 (2 9 tT {Mx)A a {V) (Mx)Y (AMT) ■ (35) 

As before, we neglect self energy corrections, so the powers of A b (x) and A c (0) must match. 

We need to solve a combinatoric problem very similar to that solved previously, namely 
the number of ways to tie the fields 


<4,m(0) pw*)T Pc(0)T> • 


(36) 
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together. We denote this as V(2n + 1,3). It can computed using our previous results for 
D{2n, 3). When n = 0, we have just a sum over one index a, so Z>(1,3) = 3. 

Now consider general n. There are two cases to consider. For the first case, A a (x) can 
tie to A a ( 0 ), 

<A,( X )A.(0)) ((AtWY (M0)T) ■ ( 37 > 

This leaves the number of ways to tie ((Ab(:r) 2 ) Tl (A c (0) 2 ) n ) together. Because of the sum 
over the index “a”, there are 3U(2n,3) ways of tying the fields together. 

The second case is if A a (x) ties to one of the other A c (0)’s. There are 2n such A c (0)’s. 
Similarly, A a (0) can then tie to one of the Ab(xy s in 2 n ways, leaving 

RW 4 ( 0 )) (A(O)AW) (A b (x)A c mA d (x) 2 r-\AM 2 ) n - 1 ) (38) 

The last correlation function is related to T>{2n — 1,3, so in total there are (2n ) 2 T>(2n — 1, 3) 
ways to tying the fields together. 

Thus we obtain 

V(2n + 1, 3) = 3 V{2n, 3) + (2 nf V(2n -1,3) =3 ( 2 n + 1)! + (2n ) 2 V(2n - 1, 3) (39) 

The solution to this equation is 

V(2n + 1,3) = (2n + 3)(2n + 1)! . (40) 

This satishes P(l, 3) = 3, and one can check directly that V( 3, 3) = 30. 

The result is then 

q 2 1 

W = (1 + z) cosh( 2 ) + (2 + z) sinh(^) , 2 = —-— . (41) 

l07T x 1 

At large distances, z —* 0 , and 

3 a 2 1 

W~l + 3z + ...~lH - — + ... , x —» cx) . (42) 

167 T xT 

This is the expected Coulomb term from single gluon exchange. 

At short distances, the exponential agrees with the result from the line-line correlator, 
Eq. (20). 

exp(z) « exp (V + log ^ + • • • , x 0 . (43) 
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The Wilson loop is similar to the exponential of the Coulomb term, but the coefficient of 
~ 1 /(xT) differs by a factor of three. It is for this reason that we refer to the behavior 
at short distances as Coulomb-like. Notice also that in the exponent, there is a factor of 
logO) r\j log (g 2 /(xT)), which has no analogy in perturbation theory. We discuss this further 
in Sec. (V). 


IV. LARGE N c 

For two colors the only difficult part of the problem was in taking into account all of 
the combinatoric factors from contracting powers of ^2 a (A a ) 2 together. For more than two 
colors, there are two difficulties. 

The first is the same problem, keeping track of combinatoric powers of ^ a (A a ) 2 , where 
the sum over the color index “a” runs from 1 to N 2 — 1. In this section we show that in 
the planar limit, taking N c —» oo while holding g 2 N c fixed, that the combinatorics is much 
simpler even than two colors. 

This simplification is special to the limit of large N c . The other problem is that for 
finite N c > 2, the symmetric structure constant, d abc , also arises in the expansion of the 
exponential. For three or more colors, there is no general expression which relates an element 
of the Lie algebra, in this case A° a) to one in the Lie group. For reasons which are clear after 
the fact, this complication can be ignored in the planar limit. 

In Appendix (A) and (B), we present the results of computations which interpolate be¬ 
tween N c — 2 and oo, for some of the lowest order coefficients which arise. They demonstrate 
the correctness of our results for N c = 2 and oo, but show that the case of general N c is not 
elementary. 

While we can compute in the planar limit, we do not find Coulomb-like behavior at short 
distances. Nevertheless, we present the results of our computation, since they are relatively 
easy to obtain. 

For three or more colors, the loop and anti-loop are not equal. We begin with the case of 
the correlator between a loop and an anti-loop, as that is more transparent. Order by order 
in perturbation theory, 

o° 1 ^ 

(L\x)L(x)) =Y. 7^5 (f) (trKW]trK( 0 )]) (44) 

n= 1 ' 
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As for two colors we neglect self energy corrections, so the number of A(rr)’s and A(0)’s must 
match. The first term, from n — 1, vanishes because of the color trace. 

For arbitary N c , the problem is computing the color traces in the above expression. Let 
us then consider the possible contractions of powers of A(x) with those of A(0). We can 
always choose the first A(x) to be left most in the trace. There are then “n” ways to connect 
this A(x ) with any of the A(0)’s. We can then make this color matrix the left most in its 
trace as well, 

V n = (tr [A n (x)} tr [A n (0)]) = n (tr [t a A n -\x)] tr [t a A n ~\0)] > . (45) 

To understand our results, consider first the case of an Abelian theory. This is much 
simpler, because then there are no color traces to bother with. For the first pair of photons, 
the factor of “n” is the same as above. For the second pair of photons, there are n — 1 ways 
to tie the second photon A{x) to one of the other n — 1 A(0)’s. This continues, so that in 
all there are n\ different ways of connecting the photons together. This cancels one of the 
factors of 1/n! in Eq. (44), giving Eq. (2). In the present case, this is 

exp ( + h wf) (46) 

This result, from single photon exchange, is valid at all distances. At large distances, this 
equals « 1 + e 2 / (AnxT) + ..., which is precisely the Coulomb term. That is, in the Abelian 
theory the Coulomb term directly exponentiates. 

In the planar limit, however, the combinatoric factors are trivial to compute. We first 
give an elementary argument, and then a detailed analysis to verify that. 

As a check, in Sec. (A) we have computed the results at finite N c for the Erst twelve 
coefficients of the perturbative expansion of the loop anti-loop correlation functions. The 
results for the first six terms agree with previous results by Burnier, Laine, and Vepsalainen 
[3] and by Brambilla, Ghiglicri, Petreczky, and Vairo [4]. Note that we only compare the 
diagrams at tree level, which are much easier than the computations to higher loop order 
performed by these authors. 

Start with Eq. (45), where we have connected one A{x) to one A(0). Consider Erst 
attaching the second gluon A(x), which is immediately to the right of t a in the color trace. 
Using the double line notation, this gives the planar diagram of Fig. (1). 

Suppose, instead, that we had attached the second gluon A{x) to the third A(0). This 
diagram is not planar, and so suppressed by 1/f/V 2 . Indeed, it is clear diagramatically that 
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FIG. 1. The planar diagram that dominates at infinite N c for the loop anti-loop correlation 
function, in the double line notation. Notice that since we use static gauge, that the ordering of 
the lines is directly in color space. 

the planar diagrams are strictly ordered in color: the second A(x) connects to the second 
.4(0), the third A(x) to the third 4(0), and so on. 

Thus in the planar limit, the combinatoric factor is simply n. This is in contrast to the 
Abelian theory, where the corresponding factor is n\, and the case of two colors, where the 
factor is (n + 1)!. 

We now proceed to an equivalent argument, which is useful in computing corrections in 
1 /N% to the planar limit. We start with the identity 

tr [fMi] [t a M 2 \ = ^tr [MiM 2 ] - trM 1 tr M 2 ■ (47) 

to reduce 

D„ = n A (tr - T tr [A"-\x)} tr p”-‘(0)]\ . (48) 

It is not evident, but one can show that the connected part of the second term is suppressed 
in the limit of large N c . 

Thus we concentrate on the first term. Note that terms of the highest order in N c originate 
from contractions as follows. Considering the term with three gluons, this is, graphically, 

A(x)A(x)A(x)A(0)A(0)A(0) (49) 
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Note that t a t a = C/I ~ N c /2 I at large N c . Using this, we work gluon by gluon to reduce 
the trace, 


(tr [A ai (x)t ai A a2 (x)t a2 A a3 (x)t a3 A bl (0)t bl A b2 {0)t b2 A b3 (0)t b3 ]) « 
A(x) (tr [t ai A a2 (x)t a2 A a3 (x)t a3 A bl (0)&A b2 (0)t b2 t ai ]) = 
A(x)C f (tr [A a 2 (x)t a 2 A a 3 (x)t a 3 A bl ( 0 )t bl A b 2 ( 0 )t b2 ]) « 

A 2 (x)Cf (tr [t°M O 3 0r)t a 3 A fel (r)Cu O2 ]) = 

A 2 (x)Cj (tr [A a 3 (x)C 3 kl 6 l (a;)t bl ]) « A 3 (x)C 3 f tr [I] = N c A 3 (x)C 3 f . 


(50) 


(51) 

(52) 

(53) 


The final factor of N c comes from the trace of a unit matrix after all contractions are done. 


In general, 


D n (x) = (tr [A n (x)] tr [A n (0)]) = — 

N c — >oo Z 



The sum of the power series gives a Bessel Function I \, 



g 2 N c l 


n 


87 r xT 


(55) 


This is our hnal result for the loop anti-loop correlator. 

Similar arguments can be used to compute the correlator between two Polyakov loops. 
This is the correlator between a (infinitely heavy) test quark and another test quark. At first 
the color structure in the planar limit looks different, but one can recognize the following. 


Take the planar diagram for the loop anti-loop correlator, and flip the line for the anti-loop 


on its head. This is then the planar diagram for the loop loop correlator. 

The only change is then the sign of the charges. For the loop anti-loop correlation 
function, each term in Eqs. (44) and (55) are positive, as the loop has charge +ig, and the 
anti-loop, charge —ig. For the loop loop correlation function, each has charge +ig, so there 
is a factor of (—l) n in the sum, 



(56) 


The propagator is a function only of the absolute value of the spatial separation, x. 
Then for the correlation functions, only the relative charges matter. Thus the correlation 
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function between a loop and an anti-loop, with charges ig and —ig, is the same as between 
an anti-loop and a loop, with charges — ig and ig, 

(■ l'Hi(O)> = (L(x)L'( 0)> . (57) 

Similarly, the correlation function between two loops, each with charge ig, equals that 
between two anti-loops, each with charge —ig, 


(L(x)L(0)) = {V(x)L-(0 )> . 


(58) 


We can now decompose the loop into real and imaginary parts. The above identities imply 
that there is no correlation between the real and imaginary parts, only between themselves. 
That for the real part is 

~l/2 

(ReL(x) ReL(O)} = 1 + ^ [/,( 2z 1 ’ 2 ) - JjpW 2 )] (59) 

and those for the imaginary parts, 

1 / 7 1 / 2 \ 

(Im L{x) ImL(O)) = ^ (^-z + — [h(2z 1/2 ) + Ji(2z 1/2 )] J (60) 

The Wilson loop can also be computed directly in the planar limit. By path ordering, 
we put all powers of A{x) to the left, and all those of .4(0) to the right. Neglecting self 
energy corrections, we expand to “n” powers of igA(x), and the same for —igA(0). There 
is a combinatorial factor of 1/n! for each term, ~ A(x) n and ~ 4(0) n . A planar diagram 
is given by connecting the left most A{x) to the left most A(0). Doing so, we then have to 
connect the second most left A{x) to the second most left 4(0). Continuing in this way, 
there is no extra factor of “n”, and the Wilson loop is given by 


1 271 1 -| / 2 7V T \ Tl 

w = nX(t) w {timx)Anm) = ^WV?r A(x) ) =Ia{2zl/2) (61) 

c n=0 y ’ n =0 v > v 7 

At large distances, z —> 0, the correlation function for the real part of the loop is 


(ReL(x)ReL(O)) « 1 + 


2iV 2 


. 2 , J 4 , , 1 (g‘N c 

z + 72 - -J “ 1 + 8ivf (‘5T 


1 

xT 


+ ... 


(62) 


The hrst term, = 1, is from the disconnected graph, which dominates at large N c . The 
second term is due to connected graphs, which are down by 1/N 2 . The coefficient of the 
term to leading order, ~ g 2 N c , agrees with the large N c limit of Eq. (19). 
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For the correlators of the imaginary part of the loop, there are no contributions from 
disconnected diagrams, only from connected diagrams, which are then necessarily r\_/ l/Nl 
For large distances, —>• 0, it begins at ~ z 3 , which is g 6 : 


(ImL(a;)ImL(0)) ~ —- ( -—h 


N} \12 2880 


g 2 N c 
96 N? V 47t 


1 

xT 


+ 


For the Wilson loop, at large distances it behaves as 
W 


z 2 £ 3 


g 2 N c 1 

1 + 2 + — + — + - — + 

4 36 87 t xT 


(63) 


(64) 


The expectation value of the Wilson loop is in and of itself a disconnected diagram, and so 
all contributions survive in the planar limit at infinite N c . The leading term, ~ g 2 , is from 
single gluon exchange. 

What is less obvious is the limit of short distances. For small z, 


Ia( z ) -t r -— 
a/27 tz 


Jl( z ) y 


sin (^) — cos ( z ) 


7 TZ 


(65) 


we hnd that the behavior of all correlation functions at small distances is not the exponential 
of ~ 1/r, but the exponential of ~ 1 j\fr. For example, for the Wilson loop 


W R 7 exp (2 a/z) rs exp 




( 66 ) 


The behavior at short distances can also be derived directly from the sum over n. At short 
distances, z is large, and one can show that the dominant behavior is given by asymptotically 
large n. To show this, one can use an effective action. For the Wilson line at large N c , it is 


where 



exp (S eff {n)) , 


S e ff(n ) = —2 log(n!) + n log z ~ n ^-oo — 2 (n log(n) — n) + n log z . 


(67) 


( 68 ) 


At large n, the stationary point of the effective action dS e ff/dn = 0 , is given by n 0 = y/z. 
Then S eff (n 0 ) 2 n 0 = 2 y/z, in accord with Eq. ( 66 ). The same method can be used 
to derive the short distance behavior of the other correlation functions. This is perfectly 
reasonable: the behavior at large distances, when z is small, is determined by the first few 
orders in perturbation theory. That at short distances, when z is large, is determined by 
the asymptotically large orders of perturbation theory. 
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V. CONCLUSIONS 


In this paper we computed the correlators of Polyakov (and Wilson) loops at tree level, 
as functions of g 2 N c / (xT). In this section we suggest that the non-Coulombic behavior at 
short distances in Eq. (66) is an artifact of the planar limit at tree level. 

We first review, on a heuristic level, why one expects Coulombic behavior at short dis¬ 
tances. As we work in imaginary time, the length in the time direction is [3 = 1/T. Now 
consider the limit of f3 —> oo. The discussion is similar for either the Wilson loop or the 
correlation function between the loop and anti-loop (or loop loop). For simplicity we discuss 
the Wilson loop, since all contributions persist in the planar limit. 

Single gluon exchange gives a contribution ~ g 2 N c /3/x. Two gluon exchange gives rise to 
two terms. The two gluons can be moved independently up and down the time axis, which 
gives the square of the first term, ~ (g 2 N c /3/x) 2 . There are also terms where the two gluons 
do not move independently, ~ ( g 2 N c ) 2 f3/x . The latter represents corrections at one loop 
order to the potential at tree level. 

For our purposes we can ignore perturbative corrections to the tree level potential. The 
leading term is from the exchange of “n” gluons, which move independently up and down the 
time axis. These one gluon kernels are identical, and so there is an associated combinatoric 
factor of 1/n!. This then generates the exponential of the potential at tree level. Because 
we assume that (3 — > oo at the outset, by dimensional analysis the coefficient of any term 
~ j3 is necessarily Coulombic, ~ 1/x. This heuristic analysis underlies extensive studies in 
perturbation theory, beginning with Ref. [7] . 

A more modern approach uses effective theories, such as potential non-relativistic QCD 
[4], For the Wilson loop, or the loop anti-loop correlator, at short distances one relates 
these correlation functions to two fields, for a color singlet and a color octet. The effective 
theory approach directly allows the exponentiation of the potential in each channel, as the 
two point function is equal to the exponential of the potential, times the expectation value 
of the associated Polyakov loop in either the singlet or octet representations. 

The same analysis can be carried out for the loop loop correlator. The operators which 
enter are those for operators with two indices in the fundamental representation. There are 
two irreducible representations, either symmetric or anti-symmetric in the two indices. (For 
example, for three colors these are the sextet and anti-triplet representations, respectively.) 
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For each representation, the two point function of the loop loop correlator is then related 
to the exponential of the associated potential, times the expectation value of the Polyakov 
loop in the associated representation. 

This analysis using effective theory applies in weak coupling for arbitrary N c , in the limit 
that f3 —> oo. This is different from our computations, where the parameter g 2 N c (3/x is 
kept finite throughout. Consequently, we believe that our result is peculiar to the planar 
approximation. We do not understand why the result at short distances is weaker than 
Coulombic in the planar limit. 

We conclude by discussing the implications of our results for the work of Lo, Friman, 
Redlich, and Sasaki [5]. Usually what is computed in non-Abelian gauge theories are the 
expectation values of (renormalized) Polyakov loops in different representations. The process 
of renormalization is not simple on the lattice [8-11]. While the self energy corrections to 
the loop vanish with dimensional regularization, they do not on the lattice. Instead, they 
give rise to terms ~ exp(— C-jzg 2 / (aT)), where “a” is the lattice spacing, and C■% the Casimir 
for the loop in an irreducible representation 1Z. By the character expansion, products of 
loops at the same point reduce to a linear sum of loops in irreducible representations [9] . 

In Refs. [5] Lo, Friman, Kaczmarek, Redlich, and Sasaki study the susceptibility of 
Polyakov loops for three colors. What is especially interesting is the difference between the 
susceptibilities for the real and imaginary parts of the loop. 

For local fields, there is no concern about the measurement of such susceptibilities. If 
one measures <p(x)(f>( 0), then by the Operator Product Expansion there are new ultraviolet 
divergences as x —> 0. However, these are local terms. A susceptibility involves the average 
of each held over the volume of space (or space-time), V. Thus while new divergences 
will arise, accompanied by powers of, e.g., the lattice spacing “a”, in the end they will be 
overwhelmed by an overall factor of 1/V. That is, the divergences associated with bringing 
two operators together does not affect the measurement of the susceptibility. 

Polyakov loops, however, are non-local operators, and on the lattice, have exponential 
divergences. Thus ultraviolet divergences, which are proportional to ~ exp(+^/(aT) m ), will 
overwhelm factors of the volume, 1/V, in the continuum limit, as aT —> 0. For operators 
which are dominated by Coulomb behavior at small distances, rn — 1. In the planar limit 
we find m = 1/2, but we believe this is special to large Ay. For the imaginary part of the 
loop, the leading perturbative correction is from three gluon exchange, Eq. (63), and so 
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then m = 3. It appears necessary to eliminate the effects of these divergences before one 
can extract susceptibilities of the Polyakov loops which are free from effects of the lattice 
cut-off. 
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Appendix A: Perturbative expansion at arbitrary N c 


In the perturbative expansion of the correlation functions, we need 


U n = <trP”(*)]taP"(0)]> = c„ (AE) 


(Al) 
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By brute force, the first eleven coefficients are 


c 2 = 


c 3 = 


c 4 = 


c 5 ~ 


C 6 = 


C 7 = 


C 8 = 

c 9 = 


N[_l 

2 2 ’ 

31V? 15 N c 

— + 

7 4 + 6 

4 

25 Nj 

+ 


2Nr 


N[ 

4 

5N[ 

32 
3 N[ 

32 

7iVj 

128 

IV? 
~~~ + 


2N? 


+ 


+ 


32 

151V? 

32 

49A? 


35N C 
8 ~ 

991V? 

16 
25971V? 


75 

4A7 


45 
+ — + 


15 

+ 7?’ 

225 225 


32 
351V 6 


128 

875A? 


+ 


8 4!V C 2 4A? ’ 

66571V, 735 2205 945 


32 
9 N? 


+ 


16 

621 Nj 


32 


+ 


64 
1545 Nj 


4 N r 161V? 


+ 


41V| 


- 903 + 


3675 2205 


21V 2 2Ah 6 


512 256 

8505 5670 


96391V 6 98731V? 


512 


128 


117451V, 90153 

' + 


16 


161V, 


21V? 


NJ 


59535 

~Wf 


(A 2 ) 

(A3) 

(A4) 

(A5) 

(A 6 ) 

(A7) 

(AS) 

(A9) 


Cio — 


Cll = 


Cl2 = 


51V 10 751V? 37051V? 94975A? 95265A? 31725 

512 32 512 256 32 4 

i 439425 212625 127575 

+ 41V? 4A? 41V? ’ 

111V ? 1 42351V? 947431V? 2051555A? 9832823A? 

2048 + 2048 + 2048 2048 + 1024 

, 13476375 4459455 12006225 8575875 779625 

+ 641V C 32A? 161V? + 161V, 7 + 4A? ’ 

31V? 2 17491V? 0 921691V? 1365573A? 2542947A? 

1024 + 1024 + 1024 1024 + 256 

, 39051045 45721665 33960465 36018675 5145525 

+ 64 161V? + 81V? + 81V? IV? 


155925 

8IV? 

(A 10 ) 


154543951V, 

256 

(AH) 


98235271V? 

128 

5145525 


41V?° 


■ (A 12 ) 


For two colors, all c n vanish for odd n, while for even n, 




(n + 1 )! 


1 

22n-2 ’ 


(A13) 


which agrees with our previous analysis: the factor of (n + 1)! is Af(n, 3) in Eq. (16). 
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Appendix B: Corrections at large N c 


From the above results, we could not deduce the general result for arbitrary n and N c . 
However, at large N c the first two terms satisfy 

Cn~n(y) (l + n ( n ~ l)K n - 1 ) -26) + 0 (J^yj ( B1 ) 

The first term, = n(N c / 2) n , agrees with our combinatorial analysis above. The second term 
is admittedly a guess, satisfied by c 2 ... c 12 . 

We now perform the exercise to assume that this guess is correct, and use it compute the 
leading correction in 1 /iV 1 2 . 

The leading order in 1/iV 2 to the correlator for the loop anti-loop correlation function is 

+ ^ 2 z2 Io(2\/z) +-^-( z - 26)Ii(2y/z)J . (B2) 

For the loop loop correlator, one changes z —> —z, 

1/7 1 z 3 4 5 6 / 2 \ 

^g ^ 2 + —£ 2 </o(2\/5) — (^ + 26) • (B3) 

Guessing the corrections to higher order is not trivial. Those ~ 1 /N^ are proportional to 
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(B4) 


when n > 4. 


[1] E. Gava and R. Jengo, Phys.Lett. B105, 285 (1981). 

[2] P. B. Arnold and L. G. Yaffe, Phys.Rev. D52, 7208 (1995), arXiv:hep-ph/9508280 [hep-ph], 

[3] Y. Burnier, M. Laine, and M. Vepsalainen, JHEP 1001, 054 (2010), arXiv:0911.3480 [hep-ph], 

[4] N. Brambilla, J. Ghiglieri, P. Petreczky, and A. Vairo, Phys.Rev. D82, 074019 (2010), 
arXiv: 1007.5172 [hep-ph], 

[5] P. M. Lo, B. Friman, O. Kaczmarek, K. Redlich, and C. Sasaki, Phys.Rev. D88, 014506 
(2013), arXiv: 1306.5094 [hep-lat]; Phys.Rev. D88, 074502 (2013), arXiv: 1307.5958 [hep-lat], 

[6] S. Borsnyi, Z. Fodor, S. D. Katz, A. Psztor, K. K. Szab, et ah, (2015), arXiv: 1501.02173 
[hep-lat] . 


22 












[7] W. Fischler, Nucl.Phys. B129, 157 (1977); T. Appelquist, M. Dine, and I. Muzinich, 
Phys.Lett. B69, 231 (1977); Phys.Rev. D17, 2074 (1978). 

[8] O. Kaczmarek, F. Karsch, P. Petreczky, and F. Zantow, Phys.Lett. B543, 41 (2002), 
arXiv:hep-lat/0207002 [hep-lat]. 

[9] A. Dumitru, Y. Hatta, J. Lenaghan, K. Orginos, and R. D. Pisarski, Phys.Rev. D70, 034511 
(2004), arXiv:hep-th/0311223 [hep-th]. 

[10] S. Gupta, K. Huebner, and O. Kaczmarek, Phys.Rev. D77, 034503 (2008), arXiv:0711.2251 
[hep-lat] . 

[11] A. Mykkanen, M. Panero, and K. Rummukainen, JHEP 1205, 069 (2012), arXiv: 1202.2762 
[hep-lat] . 


23 


